Periodic thin-film interference filters as one-dimensional photonic crystals 
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We review piiotonic band gap related properties of a simple 
periodic system of thin dielectric layers. Properties associated 
with forbidden and allowed bands of such one-dimensional 
photonic crystals are presented. A revision of forbidden bands 
properties leads to an omnidirectional Bragg mirror design. 
The anisotropy of allowed bands suggests the formation of 
photon-focusing caustics in one-dimensional photonic crys- 
tals. 



I. INTRODUCTION 

Artificial periodically microstructured dielectric mate- 
rials have attracted widespread interest in recent years 
because of the possibility to alter the dispersion rela- 
tion of photons jl]-^. Such periodic materials, pho- 
tonic crystals, can exhibit energy gaps of zero density 
of states, full photonic band gaps (PBG). After the first 
suggestions of using photonic gaps to inhibit spontaneous 
emission [^-|6| numerous attempts were made to real- 
ize three-dimensional (3D) structures that forbid visible 
light propagation in all directions (see e.g. Much 
less attention has been focused on the allowed bands 
of photonic crystals. The strong dispersion and spa- 
tial anisotropy of allowed bands lead to a number of 
new optical properties which are inconceivable in con- 
ventional crystals ]^-pl|]. In contrast to a full PBG re- 
quiring a special topology of the 3D dielectric lattice, 
as well as considerably large refractive index contrast of 
constituents, strong dispersion and anisotropy of allowed 
bands are characteristic of any periodic structures even 
with a rather small index contrast. 

The subject of this paper is the simplest type of pho- 
tonic crystals, a periodic arrangement of dielectrics layers 
of alternating refractive indeces, ni and n2, and thick- 
nesses, di and d2. Further, we use subscripts 1 and 2 for 
low and high index layers, respectively. Such kind of one- 
dimensional (ID) photonic crystals belong to the class 
of the well-known thin-film interference filters. Thin- 
film interference filters have proven to be of great im- 
portance in modern optoelectronic applications, ranging 
from Bragg mirrors for vertical cavity surface emitting 
lasers to narrow-band filters for dense wavelength divi- 
sion multiplexing systems (see e.g. |12-lJ]). However, it 



has been shown recently that some passive devices based 
on interference filters can be improved if designed on the 
basis of photonic crystals |p"5|-|l8|]. 

In general, wave propagation in periodic media can be 
described in terms of Bloch waves. For clarity, Bloch 
waves are briefly discussed in Section II. Two examples 
associated with forbidden and allowed bands of a ID 
photonic crystal are presented in Section III and Section 
IV, respectively. An example of an interference filter im- 
provement is presented in Section III. Based on photonic 
band structure analysis we show how to make a Bragg 
mirror omnidirectional. In Section IV, a discussion of 
some photon-focusing caustics patterns in ID photonic 
crystal is presented. Section V concludes the paper. 



II. BLOCH WAVES 

Optical properties of a periodic medium are described 
by its permittivity, which is a periodic function of a po- 
sition in space. In the case of ID photonic crystals a 
medium is periodic in one direction only, so the permit- 
tivity is: e{z) = e{z + /A), where z is the direction of 
periodicity, A = di -\- d2 is the period, and I is an integer. 
According to the Bloch-Floquet theorem (see e.g. [Q), 
normal electromagnetic modes of such a periodic medium 
are 



E = Eif(z) exp(— zJCz), 



(1) 
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where E/f (2) — Ea'(z + ?A) is a periodic function of pe- 
riod A. The subscript K is the Bloch wave number and 
indicates that the function E/f(z) depends on K. The 
field (^ can be expanded in Fourier series comprising 
an infinite set of partial waves, spatial harmonics. So, a 
traveling Bloch wave consists of a group of plane waves 
propagating together as a stable entity in a particular 
direction. To obtain a dispersion relation of an electro- 
magnetic wave in a periodic medium, relating the fre- 
quency of a Bloch wave, w, with the wave vector, k, a 
Fourier expansion of (|l|) should be substituted into the 
Maxwell equations. This results in a dispersion relation 
in the form of an infinite system of linear equations. The 
system should be properly truncated to get a tractable 
numerical solution. An exact analytical solution is gen- 
erally unavailable. Indeed, a ID photonic crystal is a 
unique system for which an analytical form of the disper- 
sion relation can be derived. Using the transfer matrix 
method one can obtain [||jl9|: 
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/(a;,k) = cos(ifA)- + =0. (2) 

A particular form of the functions A{u!, k^) and D{u!, k^) 
may be found elsewhere [Q^, here k^ is the tangen- 
tial component of the Bloch wave vector. The dispersion 
relation given by equation ^ contains all the informa- 
tion about photonic crystal eigenmodes, describing the 
properties of Bloch waves. 

In general, the properties of Bloch waves in a peri- 
odic structure can differ drastically from those of a plane 
waves in an isotropic homogeneous medium. Well-known 
dispersion relation of a plane wave in an isotropic homo- 
geneous nondispersive medium has a form 

f{cj,k)^{cjn/cf-k' = 0, (3) 

where n is the refractive index and c is the speed of light 
in vacuum. A 3D sketch of the dispersion relation (|^), 
usually referred to as the light cone, is depicted in the 
figure |l|(a). Only 2D slices of the wave vector space are 
presented. The main properties of the band structure 
are: (i) all polarization states are degenerate; (ii) there 
are not hand gaps; and (Hi) the direction of the energy 
flow, the direction of group velocity, Vg, points parallel to 
the wave vector k. The last property follows directly from 
the dispersion relation (|^), because the group velocity is 
the gradient of the frequency lu in the wave vector space: 

Vg = Vkuj{k). (4) 

In a periodic ID medium, even a small index perturba- 
tion yields a band structure that differs drastically from 
the light cone. In particular: (i) The polarization de- 
generacy is lifted, (ii) Photonic band gaps are devel- 
oped. In the long wavelength limit, the band structure 
asymptotically follows the light cone of some homoge- 
neous anisotropic uniaxial crystal, displaying so-called 
form birefringence For higher Bloch wave en- 

ergies, new bands, which are separated by frequency and 
angular gaps appear in the band structure [Fig. |l|(b)]. 
(Hi) The energy flow does not follow the wave vector di- 
rection any more [Fig. |^(b)]. In a periodic media, the 
velocity of the energy flow integrated over a unit cell is 
still identical to the group velocity [^. However, due 
to the strongly non-circular shape of constant-frequency 
contours in the wave vector space, the group velocity (jj) 
is no longer parallel to the wave vector. 

III. OMNIDIRECTIONAL REFLECTION 

Probably one of the most famous application of thin- 
film interference filters is a Bragg mirror, which is a peri- 
odic alternation of different dielectric layers with low and 
high indices of refraction (Fig. |2|). In general, mirrors 



come in two basic varieties: a metallic mirror with a dis- 
sipative losses of few-percent, being in fact an omnidirec- 
tional reflector and a dielectric Bragg mirror. A dielectric 
Bragg mirror can be made nearly loss-less, but it is highly 
reflecting within only limited angular range. A structure, 
combining the properties of both mirror types, i.e., being 
omnidirectional and loss-less, is of strong interest as it is 
likely to find many applications in optoelectronics and 
all-optical systems. 

Until recently, the possibility to design such a "perfect 
mirror" was mainly associated with 3D photonic crys- 
tals having a full PBG. Recently, several research groups 
worldwide have reported that a simple-to-fabricate Bragg 
mirror suffices to design a low-loss omnidirectional reflec- 
tor [p^p^ . This demonstration can lead to various high 
performance optoelectronic devices employed at any de- 
sirable wavelength. Efficient antenna substrates, energy 
saving fllters, enclosures for microcavities [^ and waveg- 
uides are a few potential applications. 

The properties of Bloch waves inside a Bragg mirror 
are governed by the dispersion relation (||). Due to the 
planar geometry of the problem, the separation of the 
electromagnetic field into TE (transverse electric) and 
TM (transverse magnetic) polarization states is possible, 
where the electric or magnetic field vector is parallel to 
the layers interfaces, respectively. This splits the prob- 
lem of light interaction with a Bragg mirror into two 
independent problems. 

A projected band structure of an infinite periodic sys- 
tem of layers is depicted in the figure ^. This is a projec- 
tion of a 3D band structure [Fig. |l|(b)] onto the lo — kx 
plane, where kx is the tangential component of the wave 
vector k, assuming that the plane of the incidence \sx — z. 
The refractive indices, ni = 1.4 and 712 = 3.4, are chosen 
close to ones of Si02 and Si in the near IR region. Thick- 
nesses of the layers are equal (c?i = o?2). The top panel 
is for TE polarization, and the bottom one for TM. An 
infinite periodic structure can support both propagating 
and evanescent Bloch waves. In figure [s], gray areas cor- 
respond to the propagating states, whereas white areas 
contain the evanescent states only and are referred to as 
photonic band gaps. 

When the frequency and the wave vector of a wave, 
impinging externally at an angle, Q;i„c, from a homoge- 
neous medium of refractive index, n, onto a thin-film 
stack, lies within the band gaps, an incident wave under- 
goes strong reflection. Pronounced high reflection bands 
(stopbands), depend strongly on frequency and angle of 
incidence. These can be easily understood from figure 
^. Photonic band gaps (i) rapidly move to higher fre- 
quencies with increasing incident angle, denoted by the 
increase of the tangential component of the wave vector; 
(ii) the TM band gaps tends to zero when approaching the 
Brewster light-line, where lu = c |kj^| /ui sinas (Fig. |^), 
as = arctan 712/711 is the Brewster angle. The TM po- 
larized wave propagates without any reflection from ni 
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to n2 layer, and from 712 to ni layer, at the Brewster an- 
gle as- These properties of the band structure restrict 
the angular aperture of a polarization insensitive range 
of high reflectance. 

In essence, omni-directional reflectance can be 
achieved due to the limitation of the number of modes 
that can be excited by externally incident waves inside 
the Bragg mirror. Light coming from the low-index am- 
bient medium (n < ^1,712) is tunneled into the internal 
cone narrowed by Snell's law (Fig. |^). Angles inside the 
crystal should be so small as to have band gaps open up 
to the grazing incident angles. In particular, (i) suffi- 
ciently large index contrast of the layers with respect to 
the ambient medium ensures that light coming from the 
outside will never go below the Brewster's angle inside 
the crystal (Fig. ^) and (ii) sufficiently large refractive 
index contrast of the layers themselves can keep the band 
gaps open up to the grazing angles |p^-p^. 

A reduced region of k-space, where electromagnetic 
modes of the photonic crystals can be excited by ex- 
ternally incident wave, lies above the light-line (Fig. ||), 
which is a 2D projection of the light cone of an ambi- 
ent medium. For example, in the case of the Si/Si02 
structure in air the first two band gaps are open for all 
external angles of incidence (shaded areas in figure ||). 
No propagating mode are allowed in the stack for any 
propagating mode in the ambient medium for either po- 
larizations. Thus, total omnidirectional reflection arises. 

To demonstrate an omnidirectional Bragg mirror, 
transmission spectra of light impinging at normal and 
oblique (30°, 60°, 89°) angles of incidence onto the 
Si/Si02 structure are shown in figure |^. Here we assume 
for calculations that ^2/(^1 + ^2) « 0.324 At nor- 
malized frequencies 0.25-0.32, a high reflectivity can be 
observed at all angles of incidence for both fundamental 
polarizations. Such a Bragg mirror can be treated as an 
omnidirectional high reflector with a relative bandwidth 
of about 25%. To obtain, e.g., an omnidirectional reflec- 
tion centered at the radiation wavelength A = 1.5/im, 
one needs a structure with a period of about 0.412/iTO. 



IV. PHOTON-FOCUSING CAUSTICS 

The strongly non-spherical shape of constant- 
frequency surfaces of ID photonic crystals [Fig. |l|(b)] 
leads to an anisotropic energy flux. We will refer to this 
phenomenon as to photon-focusing^ in analogy to phonon- 
focusing, which takes place for ballistic propagation of 
phonons in solids p^ . 

The anisotropy of periodic multilayer structures in the 
long-wavclcngth regime (A ^ A) is a well known phe- 
nomenon pO| , form birefringence. In the limit of long 
wavelengths, constant-frequency surfaces of electromag- 
netic modes in a photonic crystal are similar to the 
constant-frequency surfaces of a negative uniaxial crystal. 



At shorter wavelengths, the anisotropy becomes stronger 
and constant-frequency surfaces become distorted, es- 
pecially near the boundaries of the Brillouin zone [Fig. 
ib)]. 

The energy flow inside a photonic crystal is generally 
not along the direction of the wave vector. The energy 
propagates along the direction of the group velocity (|^), 
i.e., along the outward normal to the constant-frequency 
surface. In fact, even a small perturbation of the re- 
fractive index suffices to distort drastically a constant- 
frequency surface from its original spherical shape. As 
a consequence, an isotropic distribution of wave vectors 
emanating from an isotropic photon source inside the 
crystal does not imply an isotropic distribution of energy 
flux. 

To plot a constant-frequency surface of ID photonic 
crystals, the dispersion relation (||) should be solved at 
a fixed frequency uj. In particular, one should vary the 
value of the tangential component of the wave vector kj^ 
and solve equation (H) for a Bloch wave number K, which 
is basically the component of the Bloch wave vector. 
In figure ^ such a constant-frequency contour, which is 
an intersection of a constant-frequency surface by the 
K — kx plane, is depicted. To produce a 3D plot of a 
constant-frequency surface it is necessary to solve equa- 
tion (H) for all wave vectors in the k^ — ky plane. How- 
ever, ID photonic crystals are invariant under any rota- 
tion around the normal to the layers boundaries. That 
is why, a constant-frequency surface of a ID photonic 
crystal is simply a revolution surface of a 2D constant- 
frequency contour calculated in any plane perpendicular 
to the boundaries of the layers. 

We plot a constant-frequency contours of a homo- 
geneous isotropic (dashed curve) and a slightly modu- 
lated periodic multilayer (solid curve) media in figure 
^. The refractive index of the homogeneous medium is 
n = 3.4 and a perturbation leads to a structure with 
indeces rti — 3.4 and 712 = 3.41. We further assume 
that layers thicknesses are equal (rfi — ^2)- Constant- 
frequency contours are presented for the normalized fre- 
quency oj = 0.146, which is within a first band gap of 
the periodic medium. Near the Brillouin zone bound- 
ary the constant-frequency contour of photonic crystal 
is strongly non-circular. It consists of a concave and a 
convex regions with positive and negative Gaussian cur- 
vatures, respectively. The curvature of a photonic crys- 
tal constant-frequency contour is depicted in the inset to 
figure 1^. In contrast to the constant curvature of a homo- 
geneous medium (dashed line in the inset) the curvature 
of a photonic crystal displays a strong variation. There 
is an inficction point where the curvature vanishes (point 
A in Fig. |). 

It is known from phonon imaging experiments [ p4[ , 
that, in the geometrical optics approximation, a van- 
ishing curvature of the constant-frequency surface leads 
to an infinite energy fiux from a point source along the 
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corresponding group velocity direction. The geometri- 
cal optics approximation assumes that phonon (photon) 
wavelength is much smaller than the source and detec- 
tor sizes, and certainly much smaller than the distance 
between the source and the detector. These sharp sin- 
gularities in the energy flux are called phonon-focusing 
caustics [ p^ . 

Because the energy flux is inversely proportional to 
the curvature of the constant-frequency surfaces , the 
points with zero curvature (Fig. |^) will lead to photon- 
focusing caustics of electromagnetic energy flux inside the 
ID photonic crystal. Due to the rotation invariance of a 
photonic crystal, the Gaussian curvature of the constant- 
frequency surface vanishes along the circle produced by 
the zero curvature point A (Fig. H) . If an isotropic non- 
coherent light source is places inside a thick slab of ID 
periodic medium, then sharply deflned circular peaks in 
the intensity distribution should be detected outside the 
sample. 

To substantiate this prediction we present a Monte 
Carlo simulation of an angular intensity distribution 
which is due to an isotropic point source inside the ID 
photonic crystal. Our Monte Carlo intensity diagram 
construction process consists in generating a uniform ran- 
dom distribution of wave vectors k and solving equations 
(^ and (^) to obtain group velocities of the electromag- 
netic modes belonging to each value of k. Then the group 
velocity vectors are collected in all directions to form an 
angular plot of intensity. 

The angular distribution depicted in figure || was gen- 
erated with a large number (10^) of initial wave vectors k 
and a "source" situated inside the photonic crystal. The 
intensity was collected by a "detector" placed on the sur- 
face of the crystal along one of the layers boundary. An 
intense peak appears near the direction corresponding to 
the inflection point with zero curvature (point A in Fig. 
H). This direction corresponds to 3.6° out of the nor- 
mal to the crystal boundary. Figure ^ shows how the 
light intensity diminishes rapidly with increasing angu- 
lar deviation from that direction. At observation angles 
around zero the vanishing intensity is due to the band 
gap of the photonic crystal. 

V. CONCLUSION 

In conclusion, we have presented a brief review of the 
properties of periodic thin-film interference filters, ID 
photonic crystals. We have discussed an improvement of 
a Bragg mirror design. We have explained how to make 
a Bragg mirror omnidirectional. Finally we have pre- 
dicted the formation of photon-focusing caustics in ID 
photonic crystals. Being the simplest type of photonic 
crystals, thin-film interference filters, should be a good 
laboratory structure to study photonic band gap related 



phenomena as well as a good candidates for improve- 
ment and design of optoelectronic devices. For example, 
problems of antenna on top of an omnidirectional mirror 
and microcavity with omnidirectional Bragg mirrors are 
interesting for further study. 
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(a) Homogeneous medium 



(b) Photonic crystal 



FIG. 1. A 3D representation of the photonic band struc- 
tures of (a) an isotropic homogeneous nondispersive medium 
and (b) a ID photonic crystal. Only 2D slices of the wave 
vector space are depicted. Insets show the orientation of the 
media. A photonic crystal band structure (b) is presented 
only for one basic polarization. 




FIG. 2. Schematic representation of a dielectric multilayer 
structure. The light rays refracting and propagating through 
a stack are shown. The full domain of incident angles Oinc in 
the range from — 7r/2 to 7r/2 is mapped onto the internal cone 
of half-angle a™"^ ~ arcsinn/ni (light gray area). 




TM - polarization 



0.4 0.8 1.2 1.6 2 

Normalized wave vector 

FIG. 3. Projected band structure of a typical ID photonic 
crystal for TE (top panel) and TM (bottom panel) polariza- 
tions. The frequency and the tangential component of the 
wave vector are defined to be normalized as ujA/2nc and 
|kx|A/7r, respectively. The gray areas correspond to the 
propagating states, whereas white areas contain the evanes- 
cent states only. The shaded areas correspond to omnidi- 
rectional reflection bands. The solid lines are the ambi- 
ent-medium light-lines. 
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FIG. 6. Angular intensity distributions. Solid curve, pho- 
tonic crystal, shows an intense peake in the direction cor- 
responding to the inflection point. Dashed curve is an 
isotropic intensity distribution corresponding to a homoge- 
neous medium. 



FIG. 4. Transmission spectra from a Si/Si02 structure 
with 10 periods at different incident angles, 0°, 30°, 60°, 89°. 
The solid (dashed) curves are for TE (TM) polarization. The 
vertical lines on each spectrum are the boundaries of a spec- 
tral region when omnidirectional reflection occurs. 




FIG. 5. K — kx intersection of dispersion surfaces of homo- 
geneous isotropic (dashed curve) and slightly modulated peri- 
odic (solid curve) media. Normalized frequency is oj = 0.146. 
Layers are of equal thickness. Inset shows a local curvature 
of the dispersion contour of photonic crystal. Arrow marks 
a wave vector for which the curvature is vanishing. Concave 
and convex region of constant-frequency contours are marked 
in diflFerent shades. Vertical and horizontal axes are not to 
scale. 
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